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Spin dynamics in Heisenberg triangular antiferromagnets: A /iSR study of LiCr02 
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We report a /iSR study of LiCr02 , which has a magnetic lattice made up of a stacking of triangular 
Heisenberg antiferromagnetic (Cr"^"*", 5* = 3/2) layers. A static magnetically ordered state is observed 
below the transition temperature Tjv = 62 K, while the expected peak of the relaxation rate is 
slightly shifted downward by a few kelvins below Tjv . We draw a comparison with the isostructural 
compound NaCr02, where an exotic broad fluctuating regime has been observed [A. Olariu, P. 
Mendels, F. Bert, B. G. Ueland, P. Schiffer, R. F. Berger, and R. J. Cava, Phys. Rev. Lett. 97, 
167203 (2006)] and was suggested to originate from topological excitations of the triangular lattice. 
Replacing Na by Li strongly narrows the exotic fluctuating regime formerly observed in NaCr02, 
which we attribute to a more pronounced inter-plane coupling in LiCr02. 

PACS numbers: 76.75.+i, 74.25.Ha, 75.50.Ee 



I. INTRODUCTION 

Geometrical frustration of antiferromagnetic interac- 
tions is the object of intense research as a route to gen- 
erate novel exotic ground states. First proposed by An- 
derson for the triangular lattice spin- liquid-like be- 
havior is now commonly observed in many compounds 
with corner-shared triangular lattices. Archetypes are 
the kagome or the pyrochlore lattices 0, H, 01 • The surge 
of interest for such weakly connected lattices in the 1990s 
came from (i) the lack of good experimental realizations 
of Heisenberg triangular antiferromagnets and (ii) the 
idea that a weaker connectivity offers an additional ingre- 
dient to destabilize Neel states and generate novel exotic 
physics. 

Only very recently, the discovery of new experimental 
systems revived the interest in the pure triangular lat- 
tice. From the theoretical point of view, there was an 
early consensus that the Heisenberg triangular antiferro- 
magnet with nearest-neighbor coupling orders at T = 
in the so-called 120° configuration of spins [sj. However, 
to date, there has been little theoretical exploration in 
the dynamics associated with the triangular topology. 
Kawamura and Miyashita proposed in the 1980s that, 
even though no magnetic transition is possible at a finite 
temperature, the system goes through a cross-over regime 
due to the proliferation of topological defects Z2 . Ob- 
servations of exotic fluctuating behaviors are at the heart 
of the recent experimental explorations of the triangular 
lattice since they might have some connection with the 
aforementioned topological excitations. In NiGa2S4, a 
fluctuating behavior is observed down to low tempera- 
tures as compared to the exchange interaction 0, Q and 
recently a transition to an ordered state was finally ob- 
served at ~ 10 K The case of CS2CUCI4 [l^l is more 
controversial. Initially thought to be an anisotropic tri- 
angular lattice, the fluctuating character might be due 
to a very pronounced one-dimensional (ID) character. 



In the two previous examples, the interaction paths are 
not straightforward. NaCr02, re-discovered in the con- 
text of the surge of interest for Na cobaltates, offered a 
third, much simpler route to tackle the issue of the na- 
ture of the ground state. As shown in [ll|, it features 
well-decoupled perfect Heisenberg triangular planes. A 
broad fluctuating regime was observed well below the 
temperature at which a peak was observed in the mag- 
netic susceptibility and was tentatively associated with 
topological excitations such as Z2 vortices. 

In this paper, we present results on another stable 
member of the ACr02 family (A = Li, Na, K), namely, 
LiCr02. Our main motivation is to probe the changes 
in the exotic dynamics observed in NaCr02 through the 
modification of the structural parameters and, hence, a 
modification of the magnetic couplings. 



II. ACrOz (A=Li, Na, K) FAMILY 

As reviewed by Collins and Petrenko, most of the stud- 
ied triangular compounds had been far from ideal since 
they presented either a large anisotropy of the interac- 
tions or a sizeable single-ion anisotropy and/or next- 
neighbor coupling [T^ . To cite some of them, compounds 
from the family ABX3, where A is an alkali metal, B is a 
transition metal, and X is a halogen atom, were very well 
studied in the past. They have quasi-lD character, with 
coupling between planes much larger than the intra-plane 
one. Closer to the ideal Heisenberg case are the systems 
that belong to the family VA"2, where X=C\, Br, and I. 
However, in these compounds single-ion anisotropy was 
observed while the triangular planes are weakly coupled. 

In the ACr02 family, the magnetic structure is built by 
stacking well-separated, triangular planes of Cr"^"*" with 
spin = 3/2 in an ABC ABC sequence. A schematic 
structure is given in Fig. [T] The R3m symmetry of the 
structure ensures the isotropy of the in-plane-exchange 
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possibly antiferromagnetic in nature [l^. A ^Li NMR 
study indicated the presence, above 60 K, of a line 
shift, which demonstrates the existence of hyperfine cou- 
pling between Li and Cr, probably mediated by the oxy- 
gen [18,]. Below this temperature, a wipe out of the signal 
occurs, probably because of the fast relaxation of nuclei 
in the critical regime around the phase transition. 

In a recent NMR and /iSR study on NaCr02, an un- 
usual regime of excitations was revealed, likely related 
to the triangular nature of the magnetic structure [ll|- 
Note that no shift was found by ^"^Na measurements on 
NaCr02 in the high-T regime, which was attributed to 
the absence of hyperfine coupling between Na and Cr 
As we will detail further in this paper, the presence of 
shift on ^Li NMR spectra reveals a more substantial 3D 
coupling between layers in LiCr02. In this paper we 
present a similar /iSR study of LiCr02. Due to the 10 ns- 
15/iS window that allows detection of fast-relaxing sites, 
/iSR enables a more complete study of LiCr02 over the 
whole temperature range, than in NMR. We compare the 
temperature dependence of the relaxation rate for both 
compounds and discuss the relation with the strength of 
the coupling between planes. 

Our polycrystalline samples were synthesized by the 
stoichiometric mixture of Li2C03 and Cr203, heated at 
800°C for 48 h, with one intermediate grinding [ll|. X- 
ray measurements showed them to be single phased with 
lattice parameters a = b = 2.9 A and c — 14.4 A, in 
good agreement with previous results [l6| . 



FIG. 1: The LiCr02 structure is built of triangular planes of 
Li, O and Cr in an ABC ABC sequence. 



coupling. Also, at least in the first two members of 
the series, this coupling is restricted to the nearest- 
neighbors [ll|i 13 • Electron Spin Resonance (ESR) mea- 
surements indicate a very small single-ion anisotropy, 
thus revealing a very pronounced Heisenberg character 
which is quite common for Cr-oxides [3, fla ]. In this 
family, the c-lattice parameter increases sizably from 
14.4 A for Li to 17.9 A for K, which tunes the sepa- 
ration between Cr planes. Here the c parameter is given 
for the rhombohedral unit cell. 

An early macroscopic susceptibility study on a wide 
temperature range presented in Ref. 16] showed that all 
three compounds present antiferromagnetic correlations. 
For LiCrOa, dew = 700 K and J = 78 K, while for 
NaCrOa 0cw = 290 K and J = 40 K. 

Neutron-scattering measurements on a single crystal of 
LiCr02 [l7j show the presence of three-dimensional (3D) 
antiferromagnetic order that develops below ~ 64 K in a 
complicated 120° structure based on 18 sublattices. Ac- 
cordingly, macroscopic susceptibility and specific-heat- 
capacity data show maxima around T/v ~ 62 K [ist . Re- 
cent first-principles calculations on LiCr02 indicated the 
existence of weak coupling between the planes, which is 



III. EXPERIMENTAL RESULTS 

/xSR is a sensitive technique for the study of magnetic 
properties at a microscopic level. It relies on the implan- 
tation of positive muons, that have a spin S"^ = 1/2 and a 
very large gyromagnetic ratio, 'j/2tt = 135.5 MHz/Tesla. 
The muon spin is initially 100 % polarized opposite to 
the direction of the beam. It is commonly assumed that 
in oxides /i+ stop at ^ lA from the O^^ sites and 
the muon spin feels the dipolar magnetic field created 
by the local environment of the sample. The measured 
quantity is the asymmetry of the muon as a function of 
time, that provides information about the static and dy- 
namical magnetic behavior of the system under study. 

Complementary measurements presented in this paper 
were performed either at PSI (Switzerland) or at ISIS 
(United Kingdom) facilities. The continuous character of 
the beam at PSI enables one to get a high time resolution 
as needed at low T whereas the pulsed beam of ISIS yields 
the possibility to follow in detail slow relaxing signals as 
needed for the dynamical study. 

The typical fiSK asymmetry in zero external field (ZF), 
well below T/v, such as detected at 5 K with high statis- 
tics, is presented in Fig. [2] The spontaneous oscillations 
visible up to ~ 0.1 /zs indicate the presence of well-defined 
internal magnetic fields in the sample, showing that at 
this temperature the system is frozen and ordered. This 
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FIG. 2: Asymmetry in zero external field at early times and 
low temperature obtained for LiCr02. The spontaneous os- 
cillation is a signature of an internal field indicating ordered 
freezing at this temperature. The line is a fit to the data. In- 
set: Temperature dependance of the magnetic internal field, 
as extracted from the frequency of the oscillations. Lines are 
guides for the eye. 



finding perfectly agrees with neutron evidence [17| for a 
magnetic long-range order which can be described with 
18 sublattices of spins at 120°. In line with this quite 
involved structure, one can expect many magnetically 
inequivalent muon sites. For each site, the muon will ex- 
perience a well-defined internal field and will contribute 
to the asymmetry with a given frequency. After inspec- 
tion of the signal of Fig. [21 we were led to a minimal fit- 
ting model corresponding to a double-peaked histogram 
of the internal fields. This leads to a two-frequencies fit 
with a damping for each frequency corresponding to the 
expected but unresolved multiplicity of magnetically in- 
equivalent /i+ sites. The asymmetry was therefore fitted 
according to the formula |2C| : 



ao 



^ ^ 2 



3 z 



(1) 

Here, fi represent the weight of the two frequency com- 
ponents, found to be around 0.5 by the fit at 5 K and 
then fixed at this value at all temperatures, oq stands 
for the initial asymmetry at t = when muons are fully 
polarized (typically ao = 0.245 at PSI and oq = 0.3 at 
ISIS). The fit obtained at 5 K is represented as a con- 
tinuous line in Fig. [21 We obtained 43 and 26 MHz for 
the two frequencies of oscillations, which correspond to 
317 and 192 mT for the internal fields at the two muon 
sites. Their evolution with temperature is presented in 
the inset of Fig. [2l The field values decrease from a sat- 
uration value at low temperature down to zero around 
60 K, where the oscillation pattern becomes difficult to 
follow. The fast damping of the oscillations, taken into 
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FIG. 3: Asymmetry obtained in zero external field for the 
whole studied temperature range. Lines are for fits, as de- 
scribed in the text. 



account by the Gaussian exponential in the fitting func- 
tion, indicates a distribution of sites of width ^ 100 mT, 
in agreement with a large number of inequivalent mag- 
netic sites. 

We now turn to the study of dynamical effects, i.e., of 
electronic fluctuations, that can be tracked by the study 
of the asymmetry curves at long times. We used a com- 
bination of data from PSI and ISIS facilities. The advan- 
tage of the former is the absence of background due to 
muons not stopping into the sample, whereas the latter 
can be taken at long times with a good statistics which 
proved to be useful around the transition. Yet the un- 
certainty of the background values of the order of 0.005 
is responsible for error bars in the relaxation rate. 

Spectra obtained up to 6 fis at PSI in zero external 
field are presented in Fig. [31 In the case of a complete 
freezing, T^^ — > 0, and one obtains the so-called "one- 
third" tail at long times (a — > ao/3). When dynamical 
effects are present, fluctuations induce a relaxation of 
this tail. In the fitting formula, a stretched exponential 
with a ~ 0.5 was used empirically to account for the 
distribution of magnetically inequivalent muon sites as 
mentioned above, which yield a distribution of relaxation 
rates. 

At 5 K, the relaxation due to electronic fluctuations 
is small, Tj"^ ^ 0.03 fis^^ , which shows that the behav- 
ior is mainly static. As the temperature increases, the 
relaxation becomes more important. The one-third tail 
relaxation increases progressively up to T '--^ 56 K. At 
56 K, the tail has the fastest decrease toward zero and, 
therefore, corresponds to the maximum of relaxation. 

In the opposite limit of the fast fluctuating regime, 
above 65 K, the early time depolarization of the muon 
spin is mainly due to static nuclear dipoles, which create 
a small random fleld of ~0.4 mT, independent of temper- 
ature. The induced asymmetry has the well-established 
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Kubo-Toyabe (KT) shape, which is Gaussian-hke on the 
time scale of Fig. [31 In order to keep track of the sHght 
relaxation of electronic origin, the KT function was mul- 
tiplied by an exponential exp{-'Xt) with a weak value of 
A - 0.05 fiS-\ 

In the intermediate regime between 56 and 65 K, the 
data could not be fitted using the previous approach and 
we were led to conclude that two types of clusters, either 
static or dynamical, coexist. In order to simplify the fit- 
ting procedure, we performed experiments under a weak 
magnetic field of 3 mT, applied along the initial direction 
of the muon spin. Thus, the muon is no longer depolar- 
ized by the static field of the nuclear dipoles, which is 
much smaller than the applied field. On the other hand, 
this value is small enough so that it does not impact the 
electronic fluctuations. Asymmetry obtained at the ISIS 
facility is represented in Fig. [H We fitted the asymmetry 
for t > I according to the ansatz 



' (ao - S)ie-(*/^i)" + B, for T < 56 K 



(ao -S 



a,)ie-(*/^^) 



B + ap, 
for 56 < T < 65 K 



(ao - B)e-(*/^i)° + S, for T > 65 K 



Here the term B stands for the background due to 
muons that do not stop into the sample [B — 0.045 ± 
0.005). The term Op was introduced empirically in or- 
der to track the paramagnetic islands that survive below 
T/v, which adds to the background B. For simplicity, 
the relaxation corresponding to these sites was crudely 
neglected. The relaxation rate estimated for different 
values of ap varying from to 0.035 were used in order 
to calculate the error bars. The obtained fits are rep- 
resented as continuous lines in Fig [51 The temperature 
dependence of the relaxation rate is presented in Fig. [51 
and values obtained below 56 K are in good agreement 
with those obtained using the zero field spectra. 



IV. DISCUSSION 

The temperature dependence of the relaxation rate is 
represented in Fig. [51 and compared to the Na analog. 
The temperature scale is normalized by Tat (62 K for 
LiCr02, and 41 K for NaCr02). According to this figure, 
one can distinguish two features: 

(i) in LiCr02,the relaxation rate peak appears at ^ 
O.QTat, while in NaCr02 it appears at '-^ 0.75T/V. 

(ii) in LiCr02, the relaxation rate peak is sharper than 
for NaCr02. 

Note that in compounds that present a regular tran- 
sition, the relaxation-rate peaks where the specific heat 
does, and the width of the relaxation-rate peak is much 
smaller than in the case of NaCr02, typically 20% of 
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FIG. 4: Asymmetry obtained for an applied longitudinal field 
of 3 mT used to decouple the nuclear static contribution. 
Lines are for fits (see text). 
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FIG. 5: Relaxation rate as a function of T /Tm for NaCr02 
and LiCr02. Lines are guides for the eye. 



In NaCr02, the relaxation rate cross-over is very broad 
and appears well below T/v, which indicates the presence 
of an extended fluctuating regime. Accordingly, neutron- 
scattering results |22,] show no sign of 3D ordering above 
~ 0.75T/V but only the presence of 2D correlations. Be- 
low this temperature, an incommensurate 3D magnetic 
order develops but remains short ranged down to lowest 
temperature. 

In LiCr02, this unusual fluctuating regime appears to 
be very reduced in comparison with the Na case. The 
neutron-scattering study revealed the growing of 3D long- 
range magnetic order very close to the transition tem- 
perature [141 . Replacing Na by Li reduces the separation 
between successive Cr planes, which is likely to create a 
more important coupling in LiCr02. Note that in both 
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Na and Li compounds the dominant mechanism for in- 
plane coupUng is made by direct exchange, as shown by 
band-structure calculations [l^, as well as by calcula- 
tions that take into account the value of the Cr-O-Cr 
angle [l^]. This angle varies only little, from 96.4° in 
NaCr02 to 98.8° in LiCr02. In fact only minor changes 
occur in the Cr02 planes which cannot yield the observed 
differences between the two compounds. 

The scenario of a 3D coupling in LiCr02 is further 
supported by ''Li NMR measurements [l^, which re- 
vealed the presence of a line shift, implying the existence 
of hyperfine coupling between Li and Cr, probably me- 
diated by the oxygen. The coupling between adjacent 
Cr-^^ planes can therefore proceed indirectly by the paths 
Cr-O-Li-O-Cr or alternatively Cr-O-O-Cr. The coupling 
through direct overlapping of Cr'^+ orbitals is probably 
negligible, in view of the large separation between the 
planes, of about 4.8 A. In addition to the exchange in- 
teraction, the Cr'^^ planes are also coupled by a small 
dipolar interaction of less than 1 K. On the contrary, in 
the Na compound, the ^"^Na NMR study revealed the 
absence of line shift, and by comparison with results on 
the isostructural Naa;Co02, it was inferred that there 
is no overlapping between the Na and O orbitals plj . 
The main mechanisms for coupling, via the paths Cr-0- 
Na-O-Cr or Cr-O-O-Cr, are therefore not effective. The 
inter-layer coupling appears thus to be mainly dipolar in 
nature, so extremely small. 

The comparative study of Li and Na Chromates 
allows, thus, to track the influence of the 3D coupling on 
the magnetic properties and isolate those intrinsic to the 
triangular lattice. In NaCr02, it was suggested that the 
extended fluctuating regime is a signature of topological 
defects, Z2 vortices, whose existence was proposed long 
ago by Kawamura and Miyashita @ . The dissociation of 
vortex pairs would induce a transition at a temperature 
lower than that of the specific-heat-capacity peak. The 
present comparison underlines the sensitivity of such 
a mechanism to 3D couplings, and the sensitivity to 
similar deviations from the ideal 2D Heisenberg case 
might explain why the experimental signature of these 



defects has not been observed in real systems for a long 
time. Altogether with the recent discovery of NiGa2S4, 
this opens a new route to the exploration of the physics 
of edge-sharing triangular Heisenberg lattices. Indeed, in 
NiGa2S4, where Ni^+ ions with S — 1 form a triangular 
lattice, persisting dynamics was observed down to very 
low temperature, well below the specific-heat peak 0, 0| 
and the magnetic order revealed by neutron scattering 
measurements remains short ranged down to 0.35 K. 
Along these ideas, it would be very interesting to study 
KCr02, if this compound can be stabilized, for which 
one expects the most pronounced 2D character in 
the ACr02 family, since it has the largest separation 
between Cr planes. An early neutron-scattering work 
interestingly revealed the absence of any 3D magnetic 
ordering down to 2 K [25j . 

V. CONCLUSION 

Our //SR study shows that the triangular Heisenberg 
compound LiCr02 develops a more conventional order- 
ing than NaCr02, although the fluctuating regime re- 
mains sizeable. We propose an interpretation in terms of 
a more pronounced 3D coupling between Cr02 triangu- 
lar planes related to the smaller distance between them. 
This point is in line with neutron experiments and re- 
cent NMR results which overall reveals the importance 
of the 2D character in the study of topological excitations 
typical of triangular systems. 
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